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ABSTRACT 

The microquasar GX 339-4 experienced a new outburst in 2010: it was observed simultaneously at various wavelengths from radio up 
to soft y-rays. We focused on observations that are quasi-simultaneous with those made with the INTEGRAL and RXTE satellites: 
these were collected in 2010 March- April during our INTEGRAL Target of Opportunity program, and during some of the other 
INTEGRAL observing programs with GX 339-4 in the field-of-view. X-ray transients are extreme systems that often harbour a 
black hole, and are known to emit throughout the whole electromagnetic spectrum when in outburst. The goals of our program are 
to understand the evolution of the physical processes close to the black hole and to study the connections between the accretion 
and ejection. We analysed radio, NIR, optical, UV, X-ray and soft y-ray observations. We studied the source evolution in detail by 
producing light curves, hardness-intensity diagrams and spectra. We fitted the broadband data with phenomenological, then physical, 
models to study the emission coming from the distinct components. Based on the energy spectra, the source evolved from the canonical 
hard state to the canonical soft state. The source showed X-ray spectral variations that were correlated with changes in radio, NIR 
and optical emission. The bolometric flux increased from 0.8 to 2.9 xlO"^ erg cm"'^ s"' while the relative flux and contribution of the 
hot medium globally decreased. Reprocessing in the disc was likely to be strong at the end of our observations. The source showed a 
behaviour similar to that of previous outbursts, with some small deviations in the hard X-rays parameters' evolution. The radio, NIR 
and optical emission from jets was detected, and seen to fade as the source softened. The results are discussed within the context of 
disc and jet models. 

Key words, black hole physics - stars: individual: GX 339-4 - gamma rays: general - X-rays: binaries - infrared: general - radio 
continuum: general 

1. Introduction ing in the disc undergo inverse Comptonization. This results 

in a power law spectr um that has been detected up to a few 
X-ray transients (XTs) are accretmg low-mass X-ray binaries MeV in some sources (iGrove et aT7998': Zdziarski et al.ll2004l: 



(LMXBs) that spend most of their time in a faint, quiescent state. ICadolle Bel et alllSf iLaurent & Titarchuk .2007.1 The exact 

They undergo large amplitude outbursts with rise times of only o^jgi^ this high-energy power law spectrum is, however, sub- 

a few days or weeks (or months in the case of GX 339-4), with -^^^ ^^^^^^^ p^^j ^f^^^ high-energy emission could instead 

tynical recurrence periods of many years ( |Tanaka & Shibazaki| ^^^^^ ^^^^ synchrotron emission from a relativistic jet, that is 

[19961). It IS commonly accepted that the outburst of an XT is the usually seen i n the radio domain (iMarkoff et al.l2005l) . Note that 

result of an accretion disc instability (Frank et al. 2002). This [Laurent et al.l Eouh . through the discovery of polarized y-ray 

process triggers a transition from a quiescent state with a very e^^ssion, recently showed that both inverse Comptonisation and 

low mass accretion rate [rii « O.OIM^,/,,), where the source ac- gynchotron emission explained the 20 keV to ~ 2 MeV spectrum 

tivity is very faint or undetected (luminosity < 10^^ erg s"'), ^^e persistent X-ray binary Cygnus X-1. 

to active states. During outbursts, the optically thick and geo- ^^e spectral characteristics are coupled with different levels of 

metrically thin accretion disc has a varying inner radius (e.g., variability, quasi periodic oscillation s (QPOs) obs erved in the 

| Tomsicketal. || 2009|) and temperature, emitting at typical X-ray p^^^^. ^^^^^^^ spectrum (PDS) (e.g., iBelloni et aljl200ll l2005l) 

energies of ~1 keV. The inner regions of the disc (and the com- ^j^j^ properties of the jet emis sion a t radio, optical and 

pact object) are embedded in a hot and tenuous medium often jp, frequencies (see, e.g.. lRussell et allllo^a IC^riSet al. 20oJ). 

referred to as the corona, where soft X-ray photons originat- g^^g^ ^^e relative strengths of each spectral component, the 

degree of variability and shape of the PDS, and properties of 

Send offprint requests to: Dr. Cadolle Bel: j^e radio emission, different spectral states have been defined 
Marion.Cadolle@sciops.esa.int 
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55200 55250 55300 55350 55400 

Time (MJD) 

Fig. 1. Daily RXTE/ ASM (red), Swift/BKT (green) and INTEGRAL/lSGRl (per pointing, blue) light curves of GX 339-4 during the 
2010 outburst. Other observations performed with ATCA, the FTS, the REM/ROSS, REMIR and the ESO/ISAAC telescopes are 
indicated.The blue numbers near the bottom of the Figure are INTEGRAL revolution numbers (see Table[T]for details). 



(e.g., iMcCUntock & RemillardI 120061: iHoman & Bellonill2005h . 
The two main canonical spectral states are the low/hard state 
(LHS) and the high/soft state (HSS). In the LHS, the X-ray en- 
ergy spectrum is dominated by hard power law-like emission 
and the PDS shows a high level of variability dominated by 
a strong (~30% fractional rms) band-limited noise, associated 
with the presence of low-frequency QPOs. The LHS exhibits a 
continuously launched cor npact jet, whose radio-spectrum is flat 
(ICorbel et al.li2000. 20()1 fCiaUo et al.ll200l 120061: Ipender et all 
|2004|) . In the HSS, thermal emission from the disc dominates 
the spectrum, and the PDS shows only a weak power law noise 
component. No core radio emission is detected, and radio detec- 
tions made during the HSS are likely t o originate in discre te ejec- 
tions launched previously (see, e.g., iFender et al.ll2009l) . Other 
states have been identified and both are intermediate flavors of 
the canonical ones. These are the HIMS (hard intermediate state) 
and SIMS (soft intermediate state) as defined in Belloni (2010). 

GX 339-4 is a recurrent XT with regular outburs ts, shown to 
be a ~7 solar mass BH from d ynamical studies (iHvnes et alJ 
2003HMu noz-Darias et aD l2008h with well-sampled X-ray (e.g., 



Belloni et al. 2005) and multi-wavelength coverage during out- 



bursts (e.g.. ICorbel et^l2000l: iHoman et al.|[2005l: ICoriat et alJ 
l2009l). On Januarv 3 2010, GX 339-4 entered a new outburst 
dYamaoka eTal]|2010l) . We triggered INTEGRAL dunns the ini- 
tial hard X-ray phase (lTomsickil2010l: iPrat et all 120101) ; subse- 
quently, we triggered the second part of our INTEGRAL cam- 
paign on GX 339-4 (see Fig. [T]) in its declining hard X-ray phase 
(ICadolleBel et alJl2010h . The source increased in intensity in 



soft X-rays w hile it kept declining i n ha rd X-r ays (as announced 
in [Shaposhn ikov & Tomsickl 120101: rYull2010h . exiting the LHS 
(iBuxton et al .1120101) towards a softer state. 
Here, we report the results of our INTEGRAL obser- 
vations of GX 339-4 together with Swift, RXTE and 
UV/optical/NIR/radio data. We start with a description of the 
available data and of the analysis procedures employed in 
Sect. |2] Results are presented in Sect. [3j followed by their inter- 
pretations and modelling in Sect. H) We summarize our studies 
in Sect.|5] 



2. Observations and data reduction 

Table [1] summarizes the observations analysed in this pa- 
per, providing instrument details, energy ranges dates and 
modes. FigureHshows the RXTE/ ASM, SWIFT/BAT X-ray and 
INTEGRAL light curves of the outburst, with the time of the 
other multi-wavelength observations indicated. 



2.1. INTEGRAL 

Target of Opportunity observations were performed during 
55259-55261 and 55301-55303 MJD, corresponding to rev- 
olutions # 902 and # 916 respectively (hereafter. Rev.: an 
INTEGRAL orbit around the Earth lasts ~3 days). We also used 
data from INTEGRA L prog rams of observations of the Galactic 
Bulge ( iKuulkers et alj|2007l) . RX J1713.7-3946 and the Galactic 
disc (performed between 55237-55285 MJD, Rev. # 895-910). 



M. Cadolle Bel et al.: Multi-wavelength studies of GX 339-4 
Table 1. Log of the GX 339-4 observations analysed in this paper. 
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Observatory 


Instrument 


Bandpass 


Period 


Observation Type 








(MJD-55000) 




INTEGRAL" 


IBIS/ISGRI 


18-200 keV 


237.5-303.3 


ToO (5x5^ Hex'^) 




JEM-X 


5-25 keV 


237.5-303.3 


ToO 


RXTE" 


PCA 


3-25 keV 


208.5-345.4 


Public 


Swift'' 


XRT 


0.3-10 keV 


217.7-352.3 


Public 




BAT 


20-100 keV 


217.7-352.3 


Public 




UVOT 


180-600 nm 


217.7-352.3 


Public 


REM'' 


ROSS, REMIR 


550-2159 nm 


260.2-261.4 and 303 


ToO 




Merope 


460-883 nm 


237.7-337.8 


Public 


ATCA'' 




5.5-9 GHz 


261-263 


ToO 



" Interrupted observations: almost all Rev. (~ 3 days) between # 895-916 for INTEGRAL; Galactic Bulge monitoring for RXTE. 
* 5x5 dither pattern around the nominal target location. 
^ Hexagonal pattern around the nominal target location. 
Snap-shot observations with the specified instruments or receivers. 



The IBIS/ISGRI and JEM-X data were reduced using stan- 
dard analysis procedures of the Off-Line Scient ific Analysis 
OSA 9 .Q following standard p rocedures (e.g., Ro driguez et al.l 
120081: ICadolle Bel et all 120091) for the production of images, 
spectra and light curves. Systematic errors of 2% were added 
for both JEM-X (in the 5-25 keV range) and ISGRI (in the 18- 
400 keV range). We used the maps, the response matrices and the 
off'-axis and background corrections from OSA 9 . 0. We checked 
that the spectral index did not change by more than 2% dur- 
ing a single revolution by carrying out spectral analysis for each 
pointing. The low level of spectral variation over a single rev- 
olution allowed us to average all spectra belonging to the same 
revolution, which permitted the signal to noise ratio to be im- 
proved. The ISGRI light curves are shown in Figures [T] and |2] 
(fourth panel down). We could only perform JEM-X analysis on 
our ToO data as the source was outside the field-of-view (FOV) 
of JEM-X in the other programs. Since the RXTE/PCA instru- 
ment observed the source more frequently and had a higher sen- 
sitivity than JEM-X, the JEM-X data were not included in the 
broad-band spectra so as to be consistent over all our data sets. 
However, we verified that the best-fit spectral parameters using 
JEM-X and PCA were consistent within the error bars. 



2.2. RXTE 



2.3. Swift 

All the Swift/XRT observations discussed in this paper were 
made in Window Timing (WT) mode to avoid pile up due 
to the source brightness. Level 2 event files were produced 
following standard processing and screening procedures with 
xrtpipeline (vO.12.6). The source and background spectra 
were then extracted from these files within xselect from the 
grade events only. Source spectra were obtained within a 40- 
pixel diameter circle centred on the source position, while back- 
ground spectra were obtained using a similar sized region at an 
off-axis position. An ancillary response file taking into account 
the exposure map of the observation was generated for each 
spectrum with xrtmkarf, and the latest version of the redistri- 
bution matrix file was used in the fitting process. The original 
spectra have been rebinned by a factor of 4 before the fitting pro- 
cess. Besides, the BAT light curve (publicly available on-line) is 
shown in Figures[T]and|2](third panel down). 
In addition, data were collected from the SwiftfUVOT (UV 
and optical) instrument, with MJD spanning from 55217 to 
55310 simultaneously with INTEGRAL. The six UVOT filters 
span the wavelength range 1800-6000 A. The camera has a 
17' X 17' FOy, res ulting in images of 2048 x 2048 pixels 
(iRoming et al.l l2005h . UVOT images are combined using the 
standard UVOT uvotimsum routine and magnitudes calculated 
using the uvotsource com mand as provided by HEASARC us- 
ing an aperture of 6" (see iPoole et al.ll2008h . Detections (at a 
significance of 3cr) are discussed in Section]?] 



We analysed all available observations, taken about once ev- 
ery two days from ~ 55217 to 55311 MJD. Each observation 
lasted between ~1 and 5 ks. The RXTE data were reduced with 
the HEASO FT software package v 6.10, following standard proce- 
dures (see lRodriguez et an i2008'). Energy spectra were only ex- 
tracted from the top layer of PCA detector 2. We used the latest 
calibration files provided by the RXTE guest observer facility. A 
level of 0.6% systematic error was adopted for GX 339-4 to ac - 
count for uncertainties in the PCA response ( iJahoda et al .112006]) . 
The resultant RXTE/PCA spectra of a single observation were 
fitted simultaneously between 3-25 keV. Note that we also added 
the IBIS/ISGRI data when available; one INTEGRAL revolu- 
tion corresponds typically to three distinct RXTE observations 
(we checked the spectral parameters were not changing signifi- 
cantly and could be averaged). ASM and PCA light curves are 
shown in Figure]2](first and second panels down). 



observations were performed with 
eye moun ting ) telescope (s ee, e.g.. 



1200 U IChincarini et aP 120031: 



2.4. REMandFTS 

Optical and NIR 
the REM (r apid ey 
Zerbi & The Re m Tear^ 

Covino et al.l 120040 equipped with the ROSS optical spec- 
trograph/imager and the REMIR NIR camera. Observations 
of GX 339-4 were carried out around 55260.3, 55261.3 and 
55303.3 MJD. Image reduction was carried out by following 
standard procedures: subtraction of an averaged bias frame 
and division by a normalized flat frame. Astrometry was 
performed using the USNOBl.cQ and the 2MAS^] catalogues. 
Aperture photometry was performed with the SExtractor pack- 



http://www.nofs.navy.mil/data/fchpix/ 
http ://w w w. ipac . caltech. edu/2mass/ 
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age (iBertin & Arnoutslll996t) for all the objects in the field to 
analyse the REM data. The calibration was done against Landolt 
standard stars for the optical filters and against the 2MASS 
catalog for NIR filters. In order to minimize any systematic 
effect, we performed differential photometry with respect to a 
selection of local isolated and non-saturated standard stars. 
Optical data have also been collected using the 2-metre robotic 
Faulkes Telescope South (FTS) located at Siding Spring, New 
South Wales, Australia using the Merope camera (EM03) 
coupled with an E2V CCD42-40DD CCD. This produces a 4.'7 
X 4.' 7 FOV and images of 2048 x 2048 pixels binned 2 x 2 to 
give 1024 x 1024 pixels at 0.278 arcsec pixel"'. Science images 
were produced using an automatic pipeline which de-biases and 
flat-fields the raw imagefl Generally, 200 second exposures 
were acquired in each of the Bessel-V, Bessel-/? and SDSS-i' 
filters, approximately weekly as part of a long-ter m monitoring 
camp aign on GX 339-4 and 34 other LMXBs dLewis et al.1 
l2008h . On some dates, multiple images were taken (mostly 
in the /' band) to investigate variability on minute-timescales. 
Seeing values range from 0."8 to 3."2. Images were discarded 
if the signal-to-noise ratio was low (often due to thin cloud), or 
if the tracking or focus were poor. In total, there are 338 usable 
images taken within the time range of interest. We performed 
aperture photometry of GX 339-4 and four nearby field stars 
using APPHOT in IRAF. In addition, we used a non-variable 
star, which was fainter or of similar magnitude to GX 339-4 
to assess its level of variability. A fixed aperture radius of 6 
pixels along with a background annulus of 10-20 pixels was 
adopted for all stars in all filters. Flux calibration in V, R and ;' 
bands was achieved using photometry of the standard star field 
RU 1 49, from the list of Landolt photometric stars (iLandoltl 
which are observed regularly by FTS. We used the stars 
A, B, C, D and F as catalogued in the star field RU 149. The 
magnitudes of the stars in the RU 149 field were estimated from 
the known R and / band magnitudes using the conversion of 
iJordi et alj (l2006h . Uncertainties were calculated from the range 
of measurements from the RU 149 field stars (5 stars x 2 per fil- 
ter - 10 measurements per filter). We used the small differences 
between our measurements and these reported measurements 
to apply a small correction to the V, R and calibration of 
the GX 339-4 field stars to achieve a more accurate calibra- 
tion. FTS light curves are shown in Figure|2](lower three panels). 



2.5. Radio 

GX 339-4 was frequently observed at radio frequencies with 
the Australia Telescope Compact Array (ATCA) and the new 
CABB back-end. ATCA observed the source several times (PI 
Corbel) during the outburst. We concentrated on the simultane- 
ous observations gathered with INTEGRAL and ground-based 
instruments at 55261.89 and 55262.91 MJD. Fluxes were, at 5.5 
and 9.0 GHz, respectively, 9.02+0.10 mJy, 9.56+0.05 mJy then 
8. 16+0.05 mJy and 7.94±0. 10 mJy. 



3. Results 

The following sub-sections refer to the X/y-ray, UV, optical and 
NIR light curves shown in Figs.|2]as well as to the HID (Fig.O 
and radio data. Some spectra are shown in Figs.|4]and |5] 
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Fig. 2. RXTE/ ASM (cts s"' in 1.5-12 keV), PCA (cts s" 



http://telescope.livjm.ac.Uk/Info/TelInst/Inst/RATCam/#pipe 



in 3-30 keV), Swift/BAT (cts cm"^ s"' in 15-50 keV), 
INTEGRAL/lSGRl (cts s ' in 40-80 keV) and FTS light curves 
(filters y, R and /') obtained from 55200 to 55400 MJD. 



3.1. Light curves 

Fig. |2] presents the light curves obtained with the ASM (2-10 
keV), PCA (3-30), BAT (15-40 keV), ISGRI (40-80 keV) and 
FTS iV,R and /' magnitudes) instruments. The source shows the 
typical behaviour of a XT: the outburst started in hard X-rays 
probably around ~55199 MJD. The light curve morphologies at 
the peak of BAT and ASM are quite different. The hard X-rays 
dropped before the soft X-ray peak: around 55295 MJD, the hard 
X-rays reached their peak and the emission started to decrease 
(with a secondary peak at 55319 MJD) while the source contin- 
ued to increase in soft X-rays until 55327 MJD (i.e., more than 
a month later), before decreasing. This behaviour is indicative 
of a state transition from a hard to a softer state. The optical 
(e.g.. Fig. |2| also dropped at the start of t he transition, as re- 
ported in previou s outbursts of GX 339-4 (iHoman et al.l 120051 : 
ICoriat et al.ll2()09l) . 

While during 55259.9-55261.1 MJD (Rev. # 902) the source 
was still in a hard state, during 55301.1-55303.3 MJD (Rev. # 
916) GX 339-4 was extremely bright in the Swift/XRT (> 400 
cts s ') instrument. The JEM-X fluxes were around 370 mCrab 
and 240 mCrab in, respectively, the 3-10 and 10-20 keV energy 
bands, while the ISGRI fluxes were varying around 190 mCrab 
and 130 mCrab in the 20-40 keV and 40-80 keV energy bands. 
These numbers show the X-ray source was softening. 
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Fig. 3. HID of GX 339-4 (logarithmic scale) obtained with 
the top layer of RXTE/PCA (2.3-19.1 keV, detector 2) between 
55208 and 55347 MJDs (the numbers indicate some of the 
INTEGRAL revolutions). The source globally evolves from bot- 
tom to top, and from right to left. 
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Fig. 4. Swift/XRT (black), RXTE/PCA (red) and 
INTEGRAL/lSGRl (green) LHS count rate spectra of GX 
339-4 (first ToO of March 2010, Rev. # 902, -55259.9-55261.1 
MJD) fitted with an absorbed multi-colour disc, reflection, Fe 
line and Comptonization components. 



3.2. Hardness intensity diagram 

To study the X-ray spectral and flux evolution during the out- 
burst, we produced a hardness-intensity diagram (HID) with 
RXTE/PC A (Fig. [3l) similar t o those widely used in the litera- 
ture (e.g.. iFender et al.ll2004l) . There was no INTEGRAL cov- 
erage from the start of the outburst to the peak, hence we fo- 
cused on the main part of the outburst simultaneously with our 
INTEGRAL data (the HID extends slightly beyond the main 
data log analysed in this paper, to illustrate the evolution of the 
source after our campaign). GX 339-4 foll owed the usual pat h 
of XTs in outburst on this diagram (see, e.g.. lBelloni et al.l2005h . 
with some small deviations in the usual track. 



3.3. X-ray and y-ray spectra 

We fitted the closest (in time) XRT (0.5-9.2 keV, when avail- 
able), PC A (3-30 keV) and IBIS/ISGRI (18-200 keV) spectra 




10 

Energy (keV) 



100 



Fig. 5. Same as Fig. 31 but in distinct units: keV~ x Photons 



,-2 „-l 



keV" 



simultaneously with XSPEC vl2.6.0 ( lArnaudlll996l) . Several 
models were tested when analyzing the spectra. A normaliza- 
tion constant was added to account for uncertainties in the 
cross-calibration of the instruments. The data were in general 
well fitted using absorbed power law (or cut-off power law and 
Comptonization) with reflection combined with a multicolour 
black-body, and a Gaussian at ~ 6.5 keV attributed to a fluo- 
rescent iron (Fe) line. The line centroid could not be well con- 
strained by the PCA: we thus fixed the energy of the line to 
lie between ~ 6 keV and 7.5 keV. An additional line at 2 keV 
(due to calibration problems) was sometimes needed in XRT. 
Absorption was a co mbination of a fixed interstellar component 
(abundances fixed to IWil ms et al.l l2000l) plus a variable local 
component. XRT found the total absorption to vary between 4.9- 
6.3 X 10^' cm"^, taking into account the statistical uncertain- 
ties (+0.07), and the dependence on the position of the source 
in the HID. The lower end (4.9 x 10~' cm"^) is compatible, 
within the errors, with the average Galactic column dens ity in 
the source direction estimated from iKalberla et aP (l2005h : it is 
possible that, up to abou t 1.4 x 10^' cm"^ is intrinsic and vari- 
able (see lCabanac et al.l (l2009l) for a detailed discussion of the 
variations of the absorption in this source; XRT data were not al- 
ways available to do that in this work). For the d isc component, 
the DisKBB model in XSPEC (iMitsuda etal.lll984 was used. We 
then replaced the phenomenological models, that first allowed 
us to easily compare the spectral parameters over the outburst, 
with a more phy sical one, the thermal Comptonisation model of 
iTitarchukI (Il994l) . We tested the best modelling by adding addi- 
tional contributions and we carefully checked any improvement 
in the goodness of the fit. When a cut-off' was needed, this pro- 
vided a good fit to almost all our data (see Tableland [3] for de- 
tails; otherwise, we used a simple power law model), although 
the temperature of the seed photons (tied to the disc temperature) 
was not always well constrained. 

The best-fit parameters are reported in Table |2l Table |3] shows 
the phenomenological cut-ofF power law parameters. In addi- 
tion to providing a more physical interpretation to the data, this 
model also permits us to avoid the divergence of the power law 
flux towards low energy since Comptonisation significantly con- 
tributes above 3x kTm, and decreases significantly below this 
energy. Our spectral fitting in general lead to normalization con- 
stants very close to 1 for all instruments. Examples of the fitted 
spectra obtained during our INTEGRAL observations are shown 
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Table 2. Best-fit spectral parameters around the period of our main INTEGRAL observations. 



Time 


Observations Disc Norm.'' 


kTn, 


kT, (keV) 




iipe line 


(D''l2n 






' bol 


(MJD") 


(Rev. #) 


(keV) 


or r 




(keV) 




(dofO 





3378411™ 
35!f 
100+^g* 
100±50 

102±51 

1 74+388 

134!^:)' 
183+^f« 
370+j«^'' 
198!5f 

32816!^™. 
52!24 

840436!^ijj?^ 
174096+^j|o 

3196!^^ 
2652!|° 

QQQ+45 



217.7-217.8 
237.5-240.2 
240.7-243.2 
243.7-246.2 
246.7-249.1 
249.7-252.1 

252.5- 255.1 

255.6- 258.1 
259.9-261.1 

264.5- 267.1 

267.6- 270.1 
270.4-272.8 
273.4-276.0 
276.4-278.8 
279.4-281.5 

281.7- 281.8 
282.4-284.9 
289.3-289.6 

297.8- 297.9 
301.1-303.3 
303.6-303.8 
310.6-310.7 



N.A.s 
895'' 
896'' 
897'' 
898'' 
899'' 
900'' 
901'' 
902 
904'' 
905'' 
906'' 
907'' 
908'' 
909 
N.A.s 
910'' 
N.A.« 
N.A.s 
916 
N.A.s 
N.A.s 



0.18+0.01 
0.83+0.10 

67+" 

"•"'-0.11 

0.65±0.10 

"•" -0.07 

0.75!«i^ 
0.68t|^ 
70+" " 

"•'"-0.16 

0.19±0.01 

"•'■^-0.18 
"•"'-0.20 

0.60tli 
0.68!«?? 
62+" " 

"•"^-0.17 

21 +" "2 

"•^^-0.01 

0.85±0.09 

1 IQ+O.IO 

^•^'^-o.ii 

"•^"-0.01 

0.22±0.01 
0.61+0.02 
0.73+0.04 

^•""-0 01 



r= 1.62+0.02 

41!^ 
38+^ 
41!3 
39+^ 
40+^ 
39+^ 
45+^ 
38+3 
39+;^ 
41!^ 

38+;^ 
39+^ 
3T4 
39+3 

r= 1.80+0.01 

40!^ 
r= 1.79+0.03 
r= 1.95+0.02 

256!«3 
35^4 

p o To+0.06 



1 63+"''' 
69+81^ 

^•"^-0.14 

1 53+"i2 

I SQ+O '4 
^•-^^-0.15 
1 59+" '3 

1 59+" '3 
1 36+"'" 

^•■-^"-0.17 

1 40+0.10 

^•^^-0.11 

1 45+" " 
1 31+" '- 
1.40!°|i 
1.28t|J 
1.28t|J 
1.16t- 



•'-o.is 

I +0.14 



6.53+ 
6.56! 
6.15+ 
6.00,+ 
6.00,+"-25 
6.37+0.20 
6 17+"" 

"•-^ -0.70 

6 23+"-'' 

"•^-^-0.23 

6 23+"--' 

"•^■^-0.66 

6 27+"-" 

"•^'-0.70 

6 07+"-3" 

"•"'-0.47 

6 45+"''' 
6 29+"-35 

"•^^-0.65 

6 23+"-" 



6.27+ 
1.14+0.12 6.00,+ 



QQJ+0.02 

50+° 
"•-•"-o.so 



6 14+"-^' 
6.00,+""" 
6.00+"-22 
g-74+0.16 

6 15+^-'^^ 



34+» i» 

"•■^ -0.08 

0.25+0.05 
19+""^ 

"•^'-0.04 

0.21+0.03 
0.22+0.05 
0.26+0.05 

25+"-"5 

"•^-^-0.04 

0.32+0.04 
0.25+0.04 

34+0 05 

"•■^^-0.04 

0.35+0.04 

n Q9+0.04 
"•■^'^-0.03 

35 +" "5 
"•-^-^-0.04 

0.39+0.06 
0.39+0.04 
0.38+0.06 
0.49+0.04 
0.28+0.07 
0.10+0.05 
51+" "'' 
0.36+-« 



1.45 (281) 
0.84 (97) 
1.05 (102) 
1.13 (103) 
1.14(102) 
1.09(102) 
0.78 (102) 
0.99 (102) 

1.58 (333) 
1.07(102) 
1.11 (102) 

1.33 (102) 

1.59 (102) 
1.15 (102) 
0.88 (312) 
0.57 (255) 
0.96 (96) 

1.34 (228) 
1.32 (260) 
1.91 (252) 
1.77 (192) 
1.50(178) 



1.43 

I. 99 
2.19 
2.14 
2.29 
2.49 
2.74 
3.18 
3.38 
4.21 
4.60 
5.04 
5.67 
6.59 
7.80 
7.94 
7.69 
9.65 

II. 4 
12.2 
14.6 
13.3 



N.A. 
0.82 
0.89 
0.86 
0.92 
1.02 
1.09 
1.25 
1.64 
1.51 
1.61 
1.77 
1.98 
2.15 
2.79 
N.A. 
2.25 
N.A. 
N.A. 
2.93 
N.A. 
N.A. 



Notes: 

Models applied in XSPEC notations: constant*phabs*reflect*(diskbb-i-gaussian-i-powerlaw) or constant*phabs*reflect*(diskbb-i-gaussian-i-comptt) 
(not all components always needed). T is the photon index. A'h varied between 4.9-6.3 x 10^' cm"-. Errors are given at the 90% confidence level 
(A;r' = 2.7). 

a) MJD-55000. 

b) Disc normalization K is proportional to (R/D)- cos 6, where R is the inner disc radius in km, D is the distance to the source in kpc and 6 the 
inclination angle of the disc. 

c) Plasma optical depth. 

d) Reflection scaling factor (1 for isotropic source above disc). 

e) Computed in the 2-20 keV range. Units: 10"' erg cm"^ s"'. 

f) X-ray bands extrapolated in the 0.01 keV-10 MeV range. Units: 10" 

g) N.A. = Non Applicable (only XRT and PCA data available). 

h) XRT data not available. 

i) Pegged at lower limit (6 keV). 



' erg cm ^ s ' . 



in counts in Fig. |4] with residuals, and in Fig. |5]in energy units 
(keV^ X Photons cm"^ s"' keV"'). A small excess was seen at 
high energies, but a non-thermal modelling did not improve our 
fits: this excess was not significant. A plot summarizing the main 
spectral parameters for all the observations is shown in Fig. |6] 
Large error bars correspond to poorly constrained data. 
Apart from three very low values around 55217, 55260 and 
55280 MJDs, the disc temperature values were remarkably con- 
stant, around 0.6-0.8 keV, until the source transited to a softer 
state. After the transition, the disc temperature was first found 
at its lowest values, but it returned to its pre-transition value 
(by 55302 MJD) within a time frame of about 10-12 days, 
even though the power law index continued to increase (see 
next paragraph). The observed inner disc temperature and the 
normalization evolution which could be high (up to ~ 55303 
MJD, Rev. #916) before decreasing (see Table |2]i, correspond 
to what is usually observed for GX 339-4 in outburst (see 
Sect. 14. 111. The radius inferr ed from diskbb is however question- 
able jGierliiiski et al.ll2008l) . Several corrections need to be ap- 
plied (colour correction factor, irradiation of the disc distorting 
the spectrum, physical prescription for the inner boundary con- 
dition), but these have large uncertainties. Moreover, the size of 



the correction can vary with time. Therefore, in Fig.|6]we only 
show the fraction of the flux that is in the diskbb component 
with respect to the power law flux: as expected, its contribu- 
tion increased with time. The reflection component also evolved 
accordingly: its contribution increased with time, from 0.19 to 
0.51, compatible with the source going to a softer state. 

From 55217 to 5531 1 MJD, the power law photon index varied 
between 1.6 and 2.4 (or a cut-off power law between 1.3-2.3, 
TableO. In the early stages of the outburst, a cut-off is needed in 
the spectral fits (the Comptonization model improves the X^ed^^ 
whereas it disappears in the softer states. Indeed, later data are 
better fit with a simple power law model (Table|2]i. This indicates 
a cooling of the corona as the thermal component gets stronger, 
and then a change of the emission mechanism at the state tran- 
sition. This further confirms that the source made a tr a nsitio n 
from a hard to soft state. With the model of iTitarchukI d 19941) . 
we obtained temperatures (^T^e) and optical depth (r) values be- 
fore Rev. #916 that yield a Comptonization parameter (yxmeC^) 
of 0.10-0.22, typical for a BH in the rising phase of the LHS. 
The bolometric flux extrapolated between 0.01 keV and 10 MeV 
(note that the bulk of emission is between 0. 1 and 100 keV) rose 
from 0.8 to 2.9x 10"** erg cm"^ s"' while the 2-20 keV flux rose 



M. Cadolle Bel et al.: Multi-wavelength studies of GX 339-4 
Table 3. Cut-off power law best-fit parameters for the periods reported in Table|2] 
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Time 


INTEGRAL 


r 




2 

^red 


r,h 
P" 


(MJD") 


Observations (Rev. #) 




(keV) 


(dof) 




217.7-217.8 


N.A. 


1 S6 

-QQ3 


SOO'Ijo, 


1.56 (280) 


16 


237.5-240.2 


895 




120!}^ 
1 05+ 


0.84 (97) 


7 


240.7-243.2 


896 


1.22±0.06 


1.06 (102) 


2 


243.7-246.2 


897 


1 <-\_i_n no 
1 nn+u.uv 

^•^^-0.08 


125!^^ 


0.75 (103) 


2 


246.7-249.1 


898 


^■•'"-0.06 


1 11+'"* 


1.18 (102) 


5 


249.7-252.1 


899 


1.30±0.08 


118+1 


1.23 (102) 


3 


252.5-255.1 


900 


1.27±0.08 


108!-| 


0.77 (102 ) 


2 


255.6-258.1 


901 


1.39+0.07 


134 % 


1.02 (102) 


10 


259.9-261.1 


902 


1.69±0.04 


328!^^ 


1.69 (333) 


27 


264.5-267.1 


904 


1.30±0.04 


88±6 


1.21 (102) 


8 


267.6-270.1 


905 


1.45±0.07 




1.21 (102) 


11 


270.4-272.8 


906 


1.42±0.07 


99 


1.49 (102) 


10 


273.4-276.0 


907 


1 47+f07 




1.66(102) 


11 


276.4-278.8 


908 


1 4(^+8.8^ 

^ -^"-0.06 


gn+14 


1.17 (102) 


14 


279.4-281.5 


909 


1.77±0.04 


1 01+28 

183 2° 


0.92 (312) 


15 


281.7-281.8 


N.A. 


1.68!«;|? 


409':,4, 


1.49 (254) 


14 


282.4-284.9 


910 




126+" 


1.53 (96) 


4 


289.3-289.6 


N.A. 


1.65±0.01 


500%g 


1.73 (227) 


19 


297.8-297.9 


N.A. 


1.91±0.01 


500 %7 


1.86 (259) 


17 


301.1-303.3 


916 


2.15±0.03 


99% 


1.69 (252) 


61 


303.6-303.8 


N.A. 


1 QQ+0.08 

^^-^-0.09 


37+^0 


1.65 (192) 


84 


310.6-310.7 


N.A. 


9 ^90.06 
z..JZ._„ (17 


500 ^:,<,„ 


1.55 (177) 


100 



Notes: 

a) MJD-55000. 

b) Percentage of the disc flux versus the power law flux. 

c) Pegged at hard limit (500 keV). 



from 1.4 to 14.6 x 10 erg cm ^ s '. After that, the spectral 
parameters were compatible with the source making a transition 
from the LHS to the HIMS (until ~ 55303.6 MJD), then to softer 
states (to the SIMS until 55304.7 MJD, then quickly to the SS, 
~ 55306.1 MJD, and finally back to the SIMS and HIMS (see 
Fig.O around 55320 MJD dMottaeta l. 2010b). but this is be- 
yond the scope of our paper since it happened after the observa- 
tions listed in our Tables. 



3.4. Multi-wavelength studies 

Observations at other wavelengths of GX 339- 4 were co nducted 
with Swift/UMOT from ~ 55280 to 55300 MJD ( lYul2OI0h . aimed 
at fast optical/UV variability in the hard state and transition. No 
rapid optical/UV variability was detected due to the low UVOT 
count rates. Our long-term FTS optical light curves of this out- 
burst are shown in Fig. |2l The source brightened by ~ 4.5 mag in 
V, R and -b ands from a low lumino sity state in 2009 October, 
55126 MJD ( iLewis & Russelll l2009l) . GX 339-4 reached the 
brightest magnitudes recorded, on 2010 April 1, 55287 MJD {V 
= 15.05 ± 0.02; = 14.49 + 0.01;/' = 14.26 + 0.01), just before 
the start of the transition from the LHS to softer states. During 
this LHS rise, the source brightened at a mean rate of ~ 0.01 
mag d"' between 55237-55287 MJD. There was considerable 
optical short-term variability in the LHS. From several consecu- 
tive exposures between 55273-55320 MJD, we measured a frac- 
tional rms variability in the /' band that was fairly constant: ~ 9- 
11% with typical errors of ~ 1 % (on each date, the time reso- 
lution was between 90 and 130 seconds). This short-term vari- 
ability throughout the outburst, including the decay, will be de- 
tailed in another paper (Russell et al. 2011, in prep.). Our re- 
sults are similar to the optical fractional rms previously reported 



at ti me resolutions of < 1 sec in the LHS of GX 339-4 (see, 
e.g., iGandhi et aUlIoToh . Between 55294-55297 MJD, the op- 
tical flux faded rapidly by ~ 1 mag in 3.0 days dRussell et al.l 
I20T0I) and a change in the SED to a bluer colour was seen 
(see Sect. 14.21 ). At the same time the S'w/ff/BAT flux started to 
fade, marking the beginning of the transition away from the 
LHS (Fig. [2]). The rapid optical fading was also accompanied 
by a drop in the /'-band fractional rms variability to 1.3 + 0.7 
% by ~55301 MJD. As the source continued to soften in X-rays 
(Fig-Ell, the optical flux remained fairly constant in the following 
few months, around V ~ 16.5; R ~ 16.0 and /' ~ 16.0 mag, with 
an i'-band fractional rms of < 5 %. REM data showed that the 
source was basically at a constant luminosity on 55260-55261 
MJD (during Rev. # 902) while on 55303 MJD (during Rev. # 
916) the source faded consistently by a factor of ~ 2 in the op- 
tical (V and /?-band) and by a factor of > 5 in the NIR (/ and 
//-band), in agreement with the FTS results. Months later, af- 
ter the peak of the outburst, the optical and NIR brightnesses in 
the V, R, I, J, H and K filters all went down. A spectrum was 
taken with the ISAAC instrument on the ESO telescope between 
55261.3-55261.4 MJD and 55307.2-55307.2 MJD, and had a 
flux and spectral shape similar to our REM and FTS photomet- 
ric data (F. Rahoui et al. 201 1, in prep.). 

The ATCA observations conducted on 55283 MJD showed flux 
densities ~ 20 mJy with an inverted radio spectrum, with a spec- 
tral index in the range: -i-O.l to +0.2, typical of powerful self- 
absorbed compact jets observed in the LHS. On 55372 MJD, no 
radio emission was detected at the location of GX 339-4: this 
was consistent with the source being in a soft state (ICorbel et al.l 
i2010bi) . Ejecta very close to the core of the system were de- 
tected with a spectral index of -0.34+0.08 (which is typical of 
optically thin synchrotron emission). The new radio source was 
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Fig. 6. Summary of the parameter evolutions of GX 339-4 
along its outburst: 2-20 keV flux (10^^ erg cm"^ s"'), % of disc 
versus power law flux, high energy cut-off energy (keV), cor- 
responding photon index and inner disc temperature (keV). See 
Tables 12] and[3] and text for details. 



in the same directi on as the large sc ale jets already detected 
in GX 339-4 (see iGallo et al.l |2004|) : the impact of material 
ejected by the system into the ISM was observed (see Discussion 
Section). 



4. Discussion 

4.1. Interpretations 

According to IWu et alj (I20T0I) . this bright outburst should have 
peaked at > 0.83 Crab in hard X-rays. This is based on the 
empirical relationship between the hard X-ray peak flux and 
the waiting time since the last bright (> 0.12 Crab) outburst of 
GX 339-4. The source transited at almost the expected time and 
flux (see below); the softening of the X-ray flux at this transition 
occurred at a slightly lower flux than the 2002-2003 and 2007 
outbursts, but at a brighter flux than the 2004-2005 outburst. 
The transition may have occurred slightly earlier as there were 
a few low -level LHS mini-outbursts betw een the 2007 and 2010 
outbursts (lKongll2008l: lMarkwardlll20"09h . These mini-outbursts 
could have emptied the disc a bit more than if the source was 
in true quiescence the whole time. The BAT, ISGRI, PCA and 
ASM light curves show that the hard X-rays were dominated 
by the power law component while the soft X-rays were mainly 
due to the disc (whose flux increased with time) plus the power 
law (which was fading at the same time). When the power law 



component returns, we observe the spectral transition back from 
soft into harder states, and a secondary peak in the hard X-rays. 
The outburst evolution of GX 339-4 during the 2010 episode is 
fairly typical of BH outbursts in general. 

Our spectral analysis showed that the source transited from an 
initial LHS, with a spectrum dominated by Comptonization in 
the beginning of the outburst (~ 55217- 55281 MJD), to softer 
states where the power law cut-off was not needed anymore and 
the photon index of the power law was very soft (~ 2.3). In X- 
ray binaries, a cut-off power law spectrum is usually interpreted 
as the signature of inverse Comptonisation of soft seed photons 
by a thermalised (i.e., with velocities following a Maxwellian 
distribution) population of electrons. Changes in the hard com- 
ponent can signal the presence of a compact jet, a corona, or 
reprocessed hard X-ray emission due to X-ray heating from an 
extended central source. The hard component therefore evolved 
as expected. Regardless of the caveats mentioned in Sect. 13.31 
about the diskbb components, and the exact value of the inner 
radius obtained from the fits, we see a trend in the disc parame- 
ters (Table|2]and Fig.|6]l which is compatible with an increase of 
the mass accretion rate before the transition. Then the disc glob- 
ally recedes when th e source reaches a low luminosity (see, e.g., 
iTomsick et al.l l2009h . during the decay phase of the outburst. 
Although we have then observ ed the disc slowly moving out- 
wards during the hardening (Ch en et al.lll997tlCadolle Bel et al.l 
12004 ), this is not nece ssarily true for all transient s ources (see. 



e.g. .'Miller e t al.l2006|)- This is still strongly deb ated dDone et al.l 
2007; Rvkoff et al.ll2007l: iGierlinski et alJl2008h as, for example, 
in XTE J1817-330. We do not address this question specifically 
in this work as the majority of our data were taken in the LHS 
and not long after the main transition into softer states. In ad- 
dition, with a lower energy limit often at 3 keV, few data are 
sensitive enough for this purpose. However, one can remark that 
the rather low disc temperatures are consistent with the compact 
objec t being a BH and not a neutron star (e.g.. lTanaka & LewinI 
119951) . The soft and hard X-rays evolved as seen in previous out- 
bursts of GX 339-4: first in the HIMS, the source transited to 
softer states where the disc dominated the emission, the inner 
disc moved closer to the BH and then gradually cooled down 
before returning to quiescence. 

The spectral evolution is consistent with other works: 
RXTE monitoring showed softenin g and LFQPO ev olu- 
tions JShatJo shnikov & Tomsick 201(3 HMotta et all 1201 Oat Fv^ 
2010h as GX 339-4 started to leave the LHS jMotta et al. 

Fig . [6| sh ows some differences with the Fig. 6 of 
Motta et all (l2009l) : while the fluxes, disc temperature and 



power law indices show the same trend, our individual cut-off 
values are sometimes very high (e. g., near 55260 MJD), like in 
the HSS of their Fig. 6, and no cut-off is needed. The evolution 
in the rising phase of the LHS is therefore not as smooth as 
expected at higher energies (which are well constrained for 
the first time), but there is no clear explanation for that. The 
2-20 keV flux (from 1.43 to 14.6 x lO"** erg cm"" s"') and the 
bolometric flux (extrapolated from 0.01 to 1000 keV) variations 
(0.82 to 2.93x10"*^ erg cm"^ s"') confirm this spectral evolution. 
The source reached a high luminosity at the maximum of ~ 
12.9 X (d/6 kpc)2 10" erg s^'. This represents ~ 18% of the 
Eddington luminosity for a 7 solar mass BH; stellar mass 
BHs accreting at or below 10' ^ Lndd are found in the LHS 
dMcCUntock & Remillardll2006l) . Besides, the evolution of the 
reflection component from 0. 19 to 0.5 1 indicates that the source 
became softer Later than the observations presented in this 
paper, the source switched back to the SIMS, and subsequently 
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to the HIMS. This is not unusual: during the decay of their 
outbursts, XTE J1720-318, XTE J1650-500 and SWIFT 
J1753. 5-0127 showed a slowly receding disc with a decreasing 
inner temperature, while at the same time the relative amount 
of the power law contribution increased again. A year after the 
beginning of its outburst (end of February 2011), GX 3 39-4 
finally switched back to the LHS (Russell & Lewisir201 lah . and 
then to quiescence ( Russell & Lewis 201 Ibll 

The broadband spectral parameter evolution from hard to soft 
states, as well as our radio/N IR/optical studies, a re consistent 
with th e preliminary results of iLewis et al l (l2010h : ICorber et alJ 
(l2010ah . A discrete radio ejection usually o ccurs around the hard 
to soft transition in a BH transient outburst: ICorbel et al.l (1201 Obh 
witnessed with the ATCA the interaction of a relativistic jet from 
GX 339-4 with the interstellar medium, implying that a major 
ejection occurred earlier During the LHS, we detected the com- 
pact core jet in the o ptic al/IR. The source continued to rise in 
the LHS for a while Jwu et al.ll20Tol) . therefore it was also ris- 
ing at radio frequencies. We added radio flux data to our SED 
(Fig.|8]i but the entire campaign will be presented in a forthcom- 
ing paper (Corbel et al. 201 1, in prep.). The rapid drop in optical 
flux and colour change observed at the star t of the state transi - 
tion are reminiscent of previous outbursts dCoriat et al.ll2009h . 
Since rapid variability probably originates in the synchrotron jet 
component (Casella et al. 2010; Gandhi et al. 2010), this is con- 
sistent with the jet no longer making a st rong contribution t o 
the optical emission and that it was fading dRussell et al.ll2010l) . 
Such chromatic behaviour, also observed in the REM data, hints 
for the presence of a component mainly contributing in the NIR 
frequencies which then fades or disappears after the transition 
(see also Fig.|7]i. The optical and infrared flux continued to drop 
in the days following the transition. The changes in optical flux 
and spectru m over the transition have been repor ted before for 
GX 339-4 (iHoman et al."2005VCoriat et al.'2009) and other BH 
LMXBs (see, e.g.. iBuxton & Bailvn 2004; Russell et al. 2007,) . 
In this paper, we clearly detect the jet evolution and its dramatic 
quenching at all wavelengths. 

4.2. Spectral energy distributions 

To compute the NIR/optical spectral energy distributions (SEDs) 
of the source, we first corrected our magnitudes for interstellar 
absorption. We assum ed a colour excess of E(B -V) = 1.2 + 0.1 
dZdziarski et al.lll998h , which is fully consisten t with our mea- 
surements of A^H and of lCorbel & Fended ( 20021) . Using standard 
extinc tion curves from iFitzpatrickl d 19991) and iKataoka et al.l 
d2008l) . we obtained the de-reddening parameters for each of our 
UV/optical/NIR filters (TablelDi. Figure|7]was obtained with the 
FTS, SMARTS and REM data at four distinct epochs (see cap- 
tion). The spectra clearly changed between the FT South and 
SMARTS epochs (-55298 MJD).The jet faded over the transi- 
tion and our last REM observations (~55303 MJD), while the 
other component (probably the disc) stayed about the same. The 
i/-band faded by a factor of 10 whereas the V-band faded by a 
factor >2.5. The jet contribution moved to lower frequencies: 
data taken on ~ 55298 MJD still have the jet dominating in 
the //-band in NIR, but no longer in the optical. The last SED 
is bluer according to Fig. |7] This is exactly what one can ex- 
pect when a source undergoes a LHS to HSS spectral transi- 
tion: the colour changed and thermal processes started to domi- 
nate. During the rising LHS, Gandhi et al. (201 1, A&A, submit- 
ted) found that the jet spectrum was highly variable in the mid- 
IR. The break between optically thick (self-absorbed) and opti- 



Table 4. De-reddening parameters used to correct the 
UV/optical/NIR photometry of GX 339-4 assuming E{B-V) = 
1.21 + 0.17 mag, with the REM telescope, FTS and UVOT 



Filter 


/I, 


A,, 




(A) 


(mag) 


REM- V 


5505 


3.64 


R 


6588 


2.81 


I 


8060 


2.04 


J 


12500 


0.96 


H 


16500 


0.62 


K 


21600 


0.42 


FTS-V 


5448 


3.69 


R 


6407 


2.92 


i' 


7545 


2.27 


UV0T-MW2 


1928 


9.74 


um2 


2246 


10.97 


uwi 


2600 


7.91 


11 


3465 


5.93 


h 


4392 


4.88 


V 


5468 


3.67 



cally thin synchrotron emission in the jet spectrum was found to 
vary between ~ 3.6 and 22 fim on timescales of minutes-hours. 
Here, in Fig.|7] we are witnessing this jet component fading over 
timescales of days. 

We produced broader-band SEDs in two distinct spectral states 
(hard and soft), composed of simultaneous radio (when avail- 
able), (de-reddened) UV/optical/NIR and unabsorbed X-ray/soft 
y-ray data. They are shown in Figures [8] and |9l In Fig. [8] we plot 
the data obtained during the first INTEGRAL ToO (55259.9- 
55261.1 MJD, Rev. # 902) Fig, p is remarkably similar to the 
Fig. 2 of lCorbel & Fended (|2002|k focusing on GX 339-4 jet 
signatures during LHS. One can interpret the results with sim- 
ple power laws, but note that this is a very rough approach and 
not a physical model. An extrapolation of the radio data up to 
the NIR/optical clearly does not fit the data. Excess NIR/optical 
emission is observed. Possible sources of the residual emission 
are the disc, the irradiated face of the co mpanion star (as the 
companion canno t contribute much, see IShahbaz et al.l 1200 lb 
iHvnes et"ai] |2004|) and the jets. Similarly, several power law 
components with distinct slopes are needed to fit our SEDs. 
For example, using a simple, double power law fit to the data in 
Fig. Q gives indices of -0.9 + 0.4 for the NIR range. We fixed 
the optical component (where the disc is assumed to dominate) 
to 1.7 which is the value we found fitting a single power law 
to our SED around 55303 MJD, where the jet contribution was 
absent or still negligible (thus assuming that, at that epoch, we 
only have contribution from the disc). For this SED in the soft 
state (Fig. |9]l, we observed important spectral changes both in 
the disc and hot medium components and found that a simple 
single power law model was enough to fit the UV/optical/NIR 
data (Fig. |7]i, with an index of 1.7 + 0.2. This is indicative of 
the fading jet component in GX 339-4, as previously discussed. 
The NIR/optical ESO/ISAAC spectroscopy taken at that time 
will be commented on in a forthcoming paper (Rahoui et al., 
2011). 

After the LHS, the radio spectral index became typical of 
optically thin synchrotron radiation, probably as a result of 
freely expanding plasma blobs previously ejected (see, e.g., 
iFender et alj|2004l) . This suggests that multiple ejection events 
took place during the outburst of GX 339-4 and then interacted 
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Fig. 7. Compiled SEDs with UVOT, REM, SMARTS (M. Buxton et al. 2011, A&A, submitted) and FTS data at distinct MJD: 
MJD = 55260.21 (UVOT/mwI, u and REM/V, R, I, J, H, K); MJD = 55281.78 (UVOT/mw2 and FTS/V, R, /'); MJD = 55297.96 
(UVOT/mw2, FTS/7?, /' and SMARTS/V, H) and MJD = 55303.30 (UVOT/mw2, um2, uwl, u, b, v, SMARTS/V, /, J, H and REM/T?). 
The error bar is so large that the fluxes of uwl, um2 and uwl are comparable (within the errors). 




10-^ 10-^ 10-6 10-5 10-4 10-3 0.01 0.1 1 10 1001000 
Energy (keV) 



Fig. 8. SED from radio to soft y-ray data plotted in flux (keV^ x cm~^ x s~' x k eV~') versus energy (keV) at ~ 55259.9-55261.1 
MJD (Rev. # 902). It is remarkably similar to the Fig. 2 of lCorbel & Fendeil(l2002h . The dash-line is a simple power law extrapolation 
of the two radio data points up to the NIR frequencies (this is not a physical model, but only shown for visual purpose; see Sect. 14.2b . 



with the interstellar medium ICorbel et al.l (1201 Obi) . This could 
potentially result in re-acceleration of particles up to very 
high energies. In general, the shape of our SEDs during the 
LHS are si milar to the ones of the transient LMXB XTE 
J1118+480 jChatv et all l2003t IZurita et al.l |2006|) : in its 2000 
outburst, the SED from radio to X-rays has been explained 
as a combination of synchrotron radiation from a jet and 
a truncated optically thick disc, whereas models assuming 
advection dominated accretion flows alone underestimated 
the optical and IR fluxes dZurita et al.. ,2006t and references 



therein). In 2005, discrepancies observed between the optical 
and IR SEDs of XTE Jl 118+480 suggested that the IR was 
dominated possibly by a jet whereas the optical was dominated 
by disc emission. Power law fits to optical SEDs have also been 
performed for other BH XTs in outburst (Hynes 2005). All 
optical SEDs exhibit quasi power law spectra, with a ranging 
from 0.5-1.5, all steeper than that expected for a viscously 
heated, multi-temperature disc: Sy oc v'^^ (which differs from 
our value taken closer to the peak). Also, the authors found for 
these BHs that the UV/optical/X-ray data - when detected - 
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Fig. 9. SED from NIR to soft y-ray data plotted in flux (keV^ X cm -xs ' xkeV ') versus energy (keV) around 55301.1-55303.3 
MJD (Rev. #916). 
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Fig. 10. Unabsorbed fluxes density (in mJy) versus frequency (Hz) from the RXTE (blue), REM and UVOT (black) data of 
GX 339-4 fitted with a Gaussian, a Comptonization model plus 1) an irradiated disc (red) or 2) a viscous disc (pink) around 
55303.6-55303.8 MJD. The irradiated disc fits the data better (see Sect. 1321). 



could be fitted with a simple black-body model of an accretion the flat-spectrum synchrotron emission may be important in the 
disc heated by internal viscosity and X-ray irradiation, but the IR and optical in this source. However, they did not exclude 
inner radius could not be well constrained. They concluded that the alternative explanation that the IR excess could come from 
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the cool outer disc. More recently, iRussell et alJ (1201 ih even 
showed in the colour-magnitude diagram of XTE J 1550-564 
(Fig. 1 of their paper) that a can change dramatically over 
state transitions. In our observations of GX 339-4 presented 
here, the contribution of the radio to synchrotron emission up 
to the NlR/optical is important, and we saw it fading over the 
transition (see Sect. 14.1b . However, another component, for 
example from the cooling di sc and/or the irradiated companion , 
and/or an i rradiated disc ( Ivan Paradiis & McCUntockl [1991 
lHvnesll2005h might be necessary to account for the NIR/optical 
excess observed. Note that in quiescence, the NIR/optical of 
GX 339-4 is likely dominate d by an optically th in disc plus a 
much fainter mass donor star dShahbaz et al.ll200lT) . 



is partly based on observations with INTEGRAL, and ESA project with in- 
struments and science data centre funded by ESA member states (especially 
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(ESAC). JR acknowledges partial funding from the European Community's 
Seventh Framework Programme (FP7/2007-2013) under grant agreement num- 
ber ITN 215212 "Black Hole Universe". JAT acknowledges partial support 
from NASA under Swift Guest Observer grant NNX08AW35G. The FT North 
and South are maintained and operated by Las Cumbres Observatory Global 
Telescope Network. This research has made use of the NASA Astrophysics Data 
System Abstract Service and of the SIMBAD database, operated at the CDS, 
Strasbourg, France. 



Finally, we had sufficient NIR and optical data to ascertain when 
GX 339-4 was in the soft state and constructed the correspond- 
ing SED around 55303.7 MJD (Fig.[TOl). In fitting the SEP, we 



used a technique similar to the one described in iRahoui et al 



(I2010h: we tested a model co nsisting of a diskir (iGierliriski et al 
l2008[ iGierliilski et al l2009h added to a Gaussian, versus a model 
consisting of: diskbb + gaussian + comptt. We fixed kT^. to the 
temperature found in Table |2l We found optical slopes around 
2, which are consistent with thermal radiation. The reprocess- 
ing levels in the disc would need to be very high to explain the 
observed UV and optical fluxes. Our results are consistent with 
[Coriat et al. (2009) as they fit SEDs with a reprocessing model 
and found a peak-temperature in the UV, where the maximum 
irradiation is expected. However, we can not constrain the pa- 
rameters very well (in fact, they are underestimated) with pho- 
tometry alone. 



5. Conclusions 

We obtained a unique data set, with the largest collection 
of multi-wavelength data collected so far simultaneously for 
an outburst of GX 339-4 . We have presented light curves 
and spectral evolution of GX 339-4 in outburst from hard to 
soft states. We confirmed that the broadband behaviour was 
consistent with previous out bursts, t h ough some peculiarities 
were seen at high energies. iHvnesI (l2005h results on SEDs 
show similar behaviour f or a collection of BHs (see also, 
e.g., iRussell et aP l2006t iRussell et all l2007h . Extrapolating 
radio data to higher frequencies in the LHS SED might be 
interpreted as non-thermal emission from ejected material, 
such as discrete ejection events that rapidly faded. Jet emission 
contributed significantly to the radio, NIR and optical data. An 
optical excess was seen above the jet component and could be 
explained by contributions from the disc (the companion can not 
account much for the emission as it is fainter than R > 21 mag). 
However, alternative processes, such as X-ray irradiation, are 
sometimes inferred, depending on the source and data available. 
We now know that the amount of irradiation in NIR/optical 
varies depending on the luminosity and the X-ray state, as seen 
in our Fig. [TOl around 55303.7 MJD. We did not observe the 
source switching back to the LHS as we did fo r XTE J 1 720-3 1 8 
(ICadoUe Bel et al.lllool or XTE J1817-330 (ICadolle Bel et all 
For the first time, we have caught the jet quenching over 
tiie transition in many wavebands (Fig. |7]i. This provides an 
observational constraint on how the jet evolved over the outburst. 
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